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~\q~. Abstract 
O , 

A new mechanism for baryogenesis is proposed in the context of an extended Brans-Dicke (BD) theory. 



We generalize the BD scalar to complex field with CP violating coupling to curvature and show that the 



O 

Q • charged BD current can be enhanced during inflationary epoch. After inflation the current decays (via tree 
level interactions) into the standard model particles. When the BD scalar is charged under baryon and/or 
lepton number, the decay produces a net baryon number. Rather generically a sufficiently large scalar 
current can be produced during inflation to account for the observed baryon asymmetry of the Universe. 
Our baryogensis scenario can in an elegant way be incorporated into a model of extended inflation. 
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I. INTRODUCTION 



In an expanding universe particle densities scale rather generically as n oc 1/a 3 , where a denotes the 
scale factor of the Universe. This observation immediately implies that any charge density preceding 
the inflationary epoch is diluted during inflation into oblivion. Namely the volume during inflation 
expands so tremendously that information about any initial charge is forgotten during inflation. This 
is true of course only if the dilution is not compensated by the production. To our knowledge up to 
now no mechanism has been proposed by which a charge density in some particles would grow during 
inflation. Here we propose such a model within the context of an extended Jordan-Fierz-Brans-Dicke 
(BD) theory [lj, where the BD scalar is generalised to a complex scalar field. We shall argue that 
under rather generic conditions, an initial small scalar charge density extending over several Hubble 
volumes can be efficiently amplified during inflation. At the end of inflation the charge decays into 
ordinary matter fields, converting thus the scalar charge into baryonic or/and leptonic charges, which 
is in turn reprocessed by sphalerons into a baryon asymmetry, representing thus a mechanism for 
creation of the observed baryon asymmetry of the Universe. 

The use of scalar fields has been ubiquitous in baryogenesis models, ranging from Affleck-Dine 
baryogenesis ^S], spontaneous baryogenesis 4], to electroweak scale baryogenesis 5, 6]. Yet with a 
few exceptions [7|, |8( , inflation has not been used in connection with baryogenesis. 

The attractive feature of our baryogenesis model is that the field responsible for the strength of 
the gravitational coupling is also responsible for the generation of matter. In General Relativity no 
such mechanism is possible since the metric tensor cannot violate CP symmetry. 

The plan of the Letter is simple. In section [III we define our generalized Brans-Dicke theory. The 
resulting equations are solved and analysed in section HTT1 and in section HVl we conclude. 



II. BRANS-DICKE THEORY WITH A COMPLEX SCALAR 



In Brans-Dicke theory the gravitational constant G is promoted to a dynamical parameter changing 
in space and time l|. This variation is conventionally expressed through a scalar field $ such that 
G = Gq^ 1 with Go the bare, true constant of nature. The Brans-Dicke action is in natural units 
(c = Gq = h = 1) given by, 



S = ^ / d A x(-9) 1/2 [*R - ^(d^d^)g 



Here R denote the Ricci scalar curvature, £ ma t is the ordinary Lagrangian of matter and non- 
gravitational fields, io is the BD coupling constant and g = det((? A1!/ ). For dimensional reasons as 



well as for convenience the substitution $ — > 00 is sometimes made in literature. We will make an 
analogous substitution, since it converts the kinetic term to the canonical form. 

Here we propose the following generalisation of the Brans-Dicke theory. We promote to complex 
scalar and write the scalar sector of the theory ([1]) as follows, 

(-g)- 1/2 £t = -p(d»(j>d„W + ( - ~u;(0„^)<r + p^R + h.c.) . (2) 

Note that this lagrangian contains operators up to dimension four only, and in this sense the renormal- 
isability properties of gravity are not made worse by the theory (j2J). Furthermore, it can accomodate 
CP violation when the (dimensionless) parameters Q and p, are complex, 

Q = ue i$u , (u = \Q\), p = \p\e ie ». (3) 

There is only one physical CP- violating phase in Eq. (jSJ), however. Indeed, the phase of uj can be 
removed by a field rotation, — > 0exp(— i0 u /2), upon which the kinetic terms in Eq. (jSJ) become 
real. The remaining complex parameter, 

\x = pe~ l9uJ = fi r + ifa , (4) 

whose components are of the form, 

p r = \p\ cos(6' At - 9 W ) , m = \p\ sin(6^ - 6 J) , (5) 

mediates CP violation of the theory ([2]). Note in particular that, when 9^ — 9^ ^ rm (n =integer), then 
fii 0, and there is a nonvanishing CP violation in the theory. Moreover, when tt/2 < 9^ — 9^ < 37r/2, 
then jj, r < 0. As we shall see, the case of negative /z r 's plays an important role in the analysis that 
follows. 

Dubbing the real and imaginary component of the rotated field as <fi± such that, 

0e-tf»/2 = + + i(f)_ , (6) 

we can rewrite lagrangian (j2J) as, 

{-gY 1,2 C^ = -(p + u)(d lt <f>+d v <f>+)g^ - (p - ^{d^-d^gT + - <pl)R - 4/i,0 + 0_ J R. (7) 

The coefficients of the kinetic terms have to be such that p — uj > in order to have positive energy 
states. As can be seen the original lagrangian corresponds to two fields with different dynamics. The 
presence of charge violation should be clear by the mixing between the two fields in the last term and 
the sign difference in the quadratic terms. We emphasize that, although this model can be written 



in terms of two coupled real scalars, the physics is different. The two fields are connected through 
the relation 0e _i61 "/ 2 = <p + + icj)—, such that a non- vanishing scalar current can arise. 

Since in our model baryon number is created via a postinflationary decay of a scalar current 
j° = Q, we shall mainly be interested in the production of a scalar charge density q^ during inflation, 

g^ = ^=l(#*-0>). (8) 

A scalar field is said to be charged under the baryon number B (or any other charge not proportional 
to B+L) if - when it decays via tree level processes into the standard model matter fields - it generates 
a net lepton an baryon number according to the relation, 

Q = %V = QbB + q L L , (9) 
where Q denotes the total scalar charge, V is the spatial volume, B and L denote the baryon and 



lepton number densities, respectively, and qs and qi are the corresponding charges. Sphalerons [9 
then process any (B — L) charge produced during inflation into a net B and L according to 

8 n/ + 4(, ig + 2) 
24n / + 13(n ff + 2) V ' 

L = - 16 ^ + 9 ^ + 2 \(E-L) , (10) 
24n/ + 13(n H + 2) V 70 ' v ' 

where rip and nu denote the number of quark and lepton families and Higgs doublets, respectively. 
This can be then observed as the baryon and lepton asymmetry of the Universe today. 



III. ANALYSIS OF THE MODEL 

We shall now analyse the dynamics of the scalar field governed by the lagrangian ([2]) (and (j7j)) 
during inflation. Inflation in a Brans-Dicke theory [12] is of a powerlaw type with the Ricci curvature 
scalar decaying as R oc t~ 2 . The solutions for the scalar $ and corresponding scale factor in the 
original BD theory of Eq. ([1]) are, 

$ = M p 2 (l + X t/a) 2 
a(t) = a (l + ^/< +1/2 , (11) 

where M v = (87rG) -1 / 2 denotes the Planck mass, x 2 = S^Pz/^M 2 the Hubble constant squared in 
the Einstein frame, a 2 = (3 + 2u)(5 + 6a;), and pf is the energy density driving inflation. As can 
be seen, for early times and large u> the expansion is nearly exponential a(t) ~ e xt and R is decays 



Taction of field configurations qualify, such that the corresponding 
sj]. The amplitude of these large scale charge fluctuations, which 



slowly. For simplicity here we assume that u in Eq. (II II) is large, such that we can approximate the 
space-time by de Sitter space with the metric tensor given by 

g^ u = diag(— 1, a 2 , a 2 , a 2 ) , a = a exp(Ht) , (12) 

where a denotes the scale factor (ao is some initial scale factor of no physical relevance) and H is the 
Hubble parameter. 

In de Sitter space (T121 the lagrangian (J7|) implies the following equation of motion for the homo- 
geneous modes 0± = 4>±(t), 

4> ± + 3H<j)± t 2i ± R(j) ± + (R<j) T = , (13) 

where 

Note that we are considering the evolution of charged configurations which are to a good approxi- 
mation homogeneous at the beginning of inflation. If the Universe was highly inhomogeneous at the 
beginning of inflation, only a smal 
initial amplitude may be small 
we assume to be of statistical origin, can be estimated as follows. The number of particles N and 
antiparticles N in a volume V ~ H~ 3 are given by iV ~ T 3 V, N ~ T 3 V, and the corresponding 
charge is Q ~ iV — N, where T denotes the effective temperature (the system is not assumed to 
be in thermal equilibrium). To estimate the amplitude of charge fluctuations 8Q, note that because 
of total charge neutrality of the system, each charge is paired with an opposite charge, such that 
the boundary surface A = dV contains ~ A/£ 2 ~ T 2 V 2 ^ 3 charge pairs, where i ~ T~ l is the 
average charge separation. The statistical charge fluctuation Qa = V^a then corresponds to the 
expected net charge in the volume V, 5Q ~ Qa ~ TV 1 / 3 . Given that V 1 ^ 3 ~ l/H and T ~ ^HMp, 
we get for the charge fluctuation in a Hubble volume 5Q ~ (M P /H) 1 / 2 , which eives for the initial 
current, ~ 5Q/V ~ H 3 {Mp / H) 1 ! 2 . The charge fluctuation 8Q is of course suppressed when 
compared with the total particle number N ~ (M P /H) 3 / 2 as 5Q/N = j^/{N/V) ~ H/M P . For 
example, when this estimate is applied to the electroweak scale, one gets for the baryon to entropy 
ratio, ub/s ~ (T/Mp) 2 ~ H/M P ~ 10~ 34 . From this estimate and the analysis presented below it 
then follows that - as long as the charge does not decay during inflation - the preinflationary charge 
density fluctuations should suffice to make our baryogenesis mechanism effective. 

In de Sitter space ( fl2l) Eq. (TTBl can be solved by an exponential Ansatz. Since the equations are 
coupled, there are four independent solutions with corresponding real constants 0^]. The solutions 



read, 



3 



<t>± = & Yl ^S e ^> K ±± = ~2 H (} ± Vi-P±i) (15) 

8Rufx r 8R r~ I 17 n 

Note first that only four out of the eight constants (/>±j are independent. This is so because and 
(f>^j are related by Eq. f|T3|) as follows, 

[ K .( K . + 3 #) _ 2£ + i?]<^°j + (R^] = . (16) 

Secondly, as can easily be verified, only k and k_ + can be positive and hence can give growing 

solutions. For certain choice of parameters k can become complex, though. In order to have a 

growing current, we need two independent solutions, since the contribution of solutions of <f>± with 
equal time behavior is zero as can be seen from Eq. (jSj). The leading contribution therefore comes 

from the interference between k and Neglecting the other two decaying terms, we adopt the 

solution, 



i<g e «--t + 0ge K -+* 



(17) 



where we took the real part since <j>± are by definition real. Substituting this in definition (jSJ) of the 
current we arrive at, 

J,Hftltf!)»{(^-«--)e (K " +s - +)i } , (18) 

where we took account of the fact that in order to get a growing solution, dt[n + > 0, implying 

that must be real. 

There are two cases of interest: 

Case A. When 1 — (3 — 7 < 0, then k is complex and n_ + is real with 1 — f3 + 7 > 4, which 

implies: 

l-7</?< 7 -3, (19) 

from which we also infer that 7 > 2. In this case the current oscillates with a growing envelope 
amplitude. In order to prevent fast oscillations, some tuning is required. In particular, the 
oscillations are slow if /3 ~ 1 — 7. In this case f3 < and \i T < 0, since 7 > 2. The corresponding 
current is 



J<f> — Jq 



(i) 



cos(za_ift) H sm(ia_Ht) 



e (a +" 3)m , (20) 



where 

a + = 3 -VY^P^f, a_ = ^v^^, J « = Ha^m ~ ■ (21) 

A typical evolution of the current (|20H2Tj) is shown in figure [JJ 
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FIG. 1: The scalar charge density (current) during inflation as a function of the cosmic time (the number 
of e-foldings). Note oscillations with a growing amplitude. After inflation, the charge density freezes out. 

Case B. When 1 — (3 — 7 > 0, both k and k y are real. In this case the current grows 

provided the condition 

7 2 

P<~-J (22) 
is fulfilled. In this case no oscillations are present. From Eq. (I22p we infer that yU r < 0, and 
that there is a solution provided < 9H 2 u\/p 2 — u 2 / (4Rp). In this case the current grows 
exponentially, 

J = j(2) e (-3+a ++ a_)ff tj (23) 

with 

4 2) = - <p { :mma + - a .) . ( 24) 

A typical evolution of the current ff23l - [24|) is shown in figure [2j 

By making use of Eq. ([To]) we can relate the and (f)+l such that the prefactor in Eqs. (1211) 
and ([21]) becomes, 



J 




e-f oldings 



FIG. 2: The scalar charge density during inflation. Unlike in figure Q] no oscillations are present and the 
charge density grows exponentially. 



Note that this relation and the above analysis are valid only in the presence of CP violation when 
/ij ^ 0. Indeed, Eq. (fl3j) implies that for = the evolution of 0+ and 0_ decouple and 
and can be specified independently. In this case there is no CP violation in the lagrangian. 
Nevertheless, even in this case an initial charge can be amplified provided there is a growing solution 
for the current. The solution is however always oscillatory and oscillations are rapid. On the other 
hand, amplification within a CP violating theory is more efficient and - as we have seen above - for 
a certain choice of parameters amplification is purely exponential. 

In the above analysis we have assumed that R is varying very slowly and expansion is close to 
exponential. Most inflationary models, including inflation driven by a Brans-Dicke scalar field, are 
nearly exponential (powerlaw). Our preliminary investigation of powerlaw inflation indications that 
the current enhancement persists, albeit in that case it is weaker. 

An equivalent charge amplification can be obtained during inflation with the replacement of fi 2 R — > 
m 2 , with m being the scalar mass. In that case CP violation remains operative after inflation, 
resulting in nontrivial effects in radiation and matter epochs. On the other hand, the curvature 
coupling is turned off in postinflationary radiation era, which makes the mechanism presented here 
in some sense simpler. Moreover, in the massive scalar case, no clear connection to the gravitational 
sector exists. Thus unless one finds another motivation for the scalar field (like e.g. flat directions 
in supersymmetric affleck-Dine baryogenesis models 0, the main impetus for this scalar-field 
baryogenesis mechanism would be lost. 

After inflation the Universe enters a preheating epoch. The charged current can decay through 



a tree level coupling to ordinary matter in a charge conserving manner. Since the Brans-Dicke 
scalar is not charged under the gauge group of the Standard Model, it cannot couple directly to 
charged fermions and leptons of the Standard Model. It can however couple to neutrinos 8] via the 
Yukawa-type lagrangian, 

3 

Ant D - ipRVH^L + h - c • ( 26 ) 

Note that when <fi acquires an expectation value, this lagrangian generates Dirac neutrino masses. 
The interaction f[2"oT) induces the following tree level decay processes 

+ + (27) 

with the coupling strengths jjij and y*j, respectively, such that upon decay the scalar charge is 
converted into a lepton charge of a similar amplitude 

= Q ~ L . (28) 

When sphalerons become active (below about T ~ 10 12 GeV, which corresponds to the sphaleron 
equilibration temperature), the lepton charge (|28|) gets converted into a net baryon (and lepton) 
number according to Eq. (ITU1) . with (B — L) = —L . For example, for rif = 3 and n H = 1, Eq. (ITU1) 
gives B = -(8/37)L ~ 0.2Q. 

During preheating a lot of entropy is produced, which results in an entropy density of the order 
s ~ T 3 . Assuming instant reheating, the reheat temperature is T ~ y/HMp ~ y 1//4 , such that the 
entropy density is at most, s ~ V 3 ^ 4 . From this and the observed baryon-to-entropy ratio [3, \l4\, 

— = (8.7 ±0.3) x HT 11 , (29) 

s 

we conclude that the required scalar current at the end of inflation is of the order 

J ~ 10" V 3/4 , (30) 

where V ~ (10 16 GeV) 4 is the potential energy driving inflation. A comparison with the entropy 
density at the end of inflation, which is of the order sh ~ H 3 , is also instructive. The entropy 
production during preheating is of the order, st ~ V 3 ^ ~ (V 1 ^ / H) 3 sh ~ (Mp/H) 3 ^ 2 sh, which is 
in typical one field scalar inflationary models of the order 10 8 s#. In these models we therefore need 
to produce a scalar charge density at the end of inflation of the order, 



(31) 



When compared with the above estimate of the expected magnitude of scalar charge fluctuations in 
the primordial preinflationary chaos, the required postinflationary current (l3Tj) seems easily feasible 
within the mechanism presented here. 

One of the merits of the model is that it fits well into a model of extended inflation. Equation (jTJ) 
resembles the lagrangian of curvature coupled extended inflation as proposed by Laycock and Liddle 
except for the addition of the mixing term due to CP violation and the absence of an inflaton 
potential. Depending on the sign of /j, r , the fields <fi± assume the role of the inflaton or of the BD 
scalar. An example of the inflationary potential in which the inflaton couples to the Ricci scalar is, 

V(4>) =A(|0| 2 - |0 O |) 2 = £ A(0| - |0 O | 2 ) 2 + 2A0 2 + 2 _ (32) 

± 

where A denotes a dimensionless quartic coupling. Since the valley of the potential (1521 corresponds 
to the inflaton, for a negative fi r , 0_ corresponds to the inflaton while + is a BD scalar. As we 
argued above, this choice for \i r also provides a large parameter space for charge amplification. The 
BD scalar acquires a potential, which is however steeper than that of the inflaton. This may provide 
a way to anchoring the BD scalar to a fixed value today thus helping to make the theory compatible 
with experiments. 

IV. CONCLUSION 

By generalizing Brans-Dicke (BD) theory to a theory with a complex scalar field (T5]), we have 
constructed a baryogenesis model operative during inflation. In its most general disguise our BD 
theory contains one physical CP violating phase which can be used to amplify the corresponding 
scalar current during inflation. 

Our starting assumption is that at the beginning of inflation the Universe contains domains of 
several Hubble volumes which contain a small amount of initial scalar charge. During inflation these 
domains are stretched to a volume larger than our visible universe, and at the same time the charge 
density gets (exponentially) enhanced. When inflation ends, CP violation turnes off and the Universe 
heats up, increasing the entropy of the Universe by a large factor. The scalar charge decays in a 
charge conserving way to ordinary matter through a tree level coupling, creating a net lepton and/or 
baryon number. The lepton minus baryon number is then processed via sphalerons to the lepton and 
baryon number of a magnitude which we observe today. 

We have also shown that our mechanism can be incorporated into an extended inflationary model 
with a minimal adaptation when considering the real and imaginary part of the field as two indepen- 
dent real fields. The real part acts like a Brans-Dicke scalar while the imaginary part is the inflaton. 



The theory can provide a natural potential for the Brans-Dicke scalar, potentially anchoring the 
expectation value today and loosening the constraints on the theory. More detailed calculations are 
needed to further investigate these proposals and compare them with the experimental bounds. 

The main attractive feature of the generalised Brans-Dicke theory we consider here is in the fact 
that a model for baryogenesis is implemented in the theory of gravity itself. A related model, dubbed 
gravitational baryogenesis, was considered in Ref. [16|, where baryons are produced via a dimension 
six operator which couples a standard model fermionic current to the gradient of the Ricci scalar. The 
authors argue that in homogeneous limit this interaction can be thought of as a chemical potential 
for baryon number of the form, /xg oc R oc p 3//2 , which therefore decays as 1/a 6 in radiation era and 
as 1/a 9 / 2 in matter era. That means that, unless the Universe starts at very high energy scales in a 
very asymmetric state, that mechanism does not work. Some may consider that as a feature which is 
not very desirable. Our mechanism does not suffer from such a problem since we require only a tiny 
preinflationary charge density which extends over several Hubble volumes, and which may as well be 
generated by statistical fluctuations in preinflationary primordial chaos. 

Another interesting aspect of our model is that baryogenesis takes place during inflation, a very 
important era in the evolution of the universe not much considered in combination with baryogenesis. 
The CP and charge violation in the theory is turned off when the universe enters radiation era, such 
that after inflation the produced charge freezes in. It would be of interest to further investigate the 
model in a more general Scalar- Tensor theory of gravity or explore the analog in string theory with 
the dilaton. 
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